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In this paper we will present the self-induced approach to decoherence, which does not require 
the interaction between the system and the environment: decoherence in closed quantum systems 
is possible. This fact has relevant consequences in cosmology, where the aim is to explain the 
emergence of classicality in the universe conceived as a closed (non-interacting) quantum system. 
In particular, we will show that the self-induced approach may be used for describing the evolution 
of a closed quantum universe, whose classical behavior arises as a result of decoherence. 

> : 

I. INTRODUCTION 

' During the last years, the theory of decoherence has became the new orthodoxy in the quantum physicists commu- 
ÇS| . nity. At present, decoherence is studied and tested in many areas such as atòmic physics, quantum òptics, condensed 
matter, etc. (see references in jïj, Q). Following the initial proposal of Zeh ||, the theory was systematized and 
developed in a great number of works. According to Zurek @], ||, decoherence is a process resulting from the interac- 
tion between a quantum system and its environment; this process singles out a preferred set of states, usually called 
"pointer basis" , that determines which observables will receive definite vàlues. This means that decoherence leads to a 
sort of selection which precludes all except a small subset of the states in the Hilbert space from behaving in a classical 
7-H ' manner. Arbitrary superpositions are dismissed, and the preferred states become the candidate to classical states: 
t-H . they correspond to the definite readings of the apparatus pointer in quantum measurements, as well to the points in 
the phase space of a classical dynamical system. Environment-induced-superselection (einselection) is a consequence 
of decoherence. 

In this paper we will present a new approach to decoherence, diffcrcnt from the mainstream einselection approach. 
Oh! From this new perspective, that we will call self-induced approach, decoherence does not require the interaction 
+1> ■ between the system and the environment: decoherence in closed quantum systems is possible. This fact has relevant 
consequences in cosmology, where the aim is to explain the emergence of classicality in the universe conceived as a 
closed (non-interacting) quantum system. In particular, we will show that the self-induced approach may be used for 
describing the evolution of a closed quantum universe, whose classical behavior arises as a result of decoherence. 
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II. THE EINSELECTION APPROACH IN COSMOLOGY 

According to Zurek, einselected states are distinguished by their stability in spite of the monitoring environment. 
In Paz and Zurek's words M, "the environment distills the classical essence of a quantum system". This 
means that, from the einselection view, the split of the universe into the degrees of freedom which are of direct 
interest to the observer (the system) and the remaining degrees of freedom (the environment) is absolutely essential 
for decoherence. Such a split is necessary, not only for explaining quantum measurement, but also for understanding 
que quantum origin of the classical world. In fact, the einselection approach considers the problem of the transition 
from quantum to classical as the core of the problem: quantum measurement is conceived as a particular case of 
the general phcnomenon of the emergence of classicality. In addition, if classicality only emerges in open quantum 
systems, it must always be accompanied by other manifestations of openness, such as dissipation of energy into the 
environment. Zurek [[2| even considers that the prejudice which seriously delayed the solution of the problem of the 
transition from quantum to classical is itself rooted in the fact that the role of the openness of a quantum system in 
the emergence of classicality was ignored for a very long time. 

In summary, decoherence explains the emergence of classicality, but only open systems can decohere. The question 
is: what about the universe as a wholc? If, as Zurek himsclf admits, the universe is, by definition, a closed system, then 
it cannot decohere. How to explain, then, the classical behavior of stars, galàxies and clusters? Zurek Q considers 
this possible criticism: if the universe as a whole is a single entity with no "outside" environment, any resolution 
involving its division into systems seems unacceptable. Zurek's answer to this objection is based on his particular 
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conception about the nature of quantum mechanics: for him, the aim of the theory is to establish the relationships 
between formal results and observer's perceptions. And perception is an information processing function carried out 
by a physical system, the brain. The brain is conceived as a massive, neural network-like computer very strongly 
coupled to its environment, and the environment plays the role of a commonly accessible intcrnet-likc data base, which 
allows the observer to make copies of the records concerning the states of the system with no danger of altering it [p] . 
The stability of the correlations between the state of the observer's brain and the state of the environment on the one 
hand, and between the state of the environment and the state of the observed system on the othcr, is rcsponsible for 
the perception of classicality. 

This means that, for the einselection approach, the problem of the transition from quantum to classical amounts to 
the question "why we don't perceive superpositions?" ||. In other words, the task is to explain, not the emergence of 
classicality, but our perception of classicality. But this position would hardly convince the cosmologist, who conceives 
the universe as a single closed object with no other object to interact with. In the cosmological context, the wave 
function of the universe describes, not the system of everything except the observers' brains, but the universe as a 
whole. Nevertheless, cosmology tries to explain, with the universal wave function, the evolution of a closed quantum 
universe whcrc the classical behavior described by general relativity emerges. If we take Zurek's position seriously, 
without the assumption of a division of the universe into individual systems the problem of the emergence of classicality 
has no solution. 

At this point, it could be noted that the einselection approach has been applied to the cosmological level with 
interesting results. This is certainly true, but does no undermine the closed-universe objection. In the works whcre 
the einselection approach is used in cosmology, the general strategy consists in splitting the universe into some degrees 
of freedom which represent the "system" of interest, and the remaining degrees of freedom that are supposed to be 
non accessible and, therefore, play the role of an internal environment. For instance, in quantum field theory, it is 
usual to perform a decomposition on a scalar field tf>, <j> = <f>s + where <j)s denotes the system field and <pE denotes 
the environment field; when it is known that the background field follows a simple classical behavior, the scalar field is 
decomposed according to <j> — (f) c + <f> q , where the background field (f) c plays the role of the system and the fluctuation 
field 4>q plays the role of the environment (see 0). This means that, strictly speaking, it is not the universe what 
decoheres, but a subsystem of the universe: we perceive a classical universe because there are unaccessible degrees of 
freedom that act as an environment. 

These considerations allows us to point out the weakest spot of the einselection program. When this approach 
is applied to the universe (and, in general, to any system with internal environment), the space of the observables 
which will behave classically is assumed in advance: the distinction between the system's degrees of freedom and 
the environmental degrees of freedom is established in such a way that the system decoheres in some observables of 
that space. This means that the split of the whole must be decided case by case: there is not a general criterion 
for discriminating between system and environment. In fact, in the case of the decomposition of the scalar field <p 
previously mentioned, different criteria are used: sometimes the decomposition is performed on the basis of the length, 
mass or momentum scales of the system and the environment, sometimes the system field is considered as containing 
the lower modes of <p and the environment as containing the higher modes (see 0). Zurek |6) recognizes that this lack 
of a general criterion for deciding where to place the "cut" between system and environment is a serious difficulty 
of his proposal: "In particular, one issue which has been often taken for granted is looming big, as a 
foundation of the whole decoherence program. It is the question of what are the 'systems' which play 
such a crucial role in all the discussions of the emergent classicality. This issue was raised earlier, but 
the progress to date has been slow at best". 

As we will see, the self-induced approach to decoherence overcomes these problems to the extent that it does not 
require the openness of the system of interest and its interaction with the environment. 

III. THE SELF-INDUCED APPROACH TO DECOHERENCE 

This approach relies on the general idea that the interplay between observables and states is a fundamental element 
of quantum mechanics (see [|[). The departing point consists in adopting an àlgebra of observables A as the primitive 
element of the theory: quantum states are represented by linear functionals over A. In the original formulation of the 
algebraic formalism, the àlgebra of observables is a C*-algebra. The GNS theorem (Gel'fand-Naimark-Segal) proves 
that the traditional Hilbert space formalism is a particular representation of this algebraic formalism; the àlgebra of 
observables is thereby given a concrete representation as a set of self-adjoint bounded operators on a separable Hilbert 
space. Nevertheless, it is well known that the C*-algebraic framework does not admit unbounded operators; therefore, 
it is necessary to move to a less restrictive framework in order to accommodate this kind of operators. The self-induced 
approach adopts a nuclear àlgebra B as the àlgebra of observables A: its elements are nuclei or kernels, that is, two 
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variables distributions triat can be though of as generalized matrices |Ï(J . By means of a generalized version of the 
GNS theorem jïïj it can be proved that this nuclear formalism has a representation in a rigged Hilbert space: the 
appropriate rigging provides a mathematical rigorous foundation to unbounded operators (see 0]). In fact, the 
nuclear spectral theorem of Gel'fand and Maurin establishes that, under very general mathematical hypothcses (quite 
reasonable from a physical point of view), for every CSCO (completo set of commuting observables) of essentially self- 
adjoint unbounded operators, there is a rigged Hilbert space where such a CSCO can be given a generalized eigenvalue 
decomposition, meaning that a continuum of generalized eigenvalues and eigenvectors may thereby be associated with 
it. In order to find the appropriate rigging, the nuclear àlgebra is used to generate two additional topologies: one of 
them corresponds to a nuclear space, which is the space of generalized observables Vb; the other corresponds to the 
dual of the space Vo, and this is the space Vs of states. 

Following ||], |Ï3| ], [H, we will symbolize an observable belonging to Vo by a round ket \0) and a state belonging 
to Vs by a round bra [p\. The result of the action of the round bra (p\ on the round ket \0) is the mean value of the 
observable \0) in the state (p\ : 

(0) P = ( P \0) (1) 

If the basis is discrete, (0) p can be computed as usual, that is, as Tr(pO). But if the basis is continuous, Tr(pO) is 
not well defined; nevertheless, (p\0) can always be rigorously defined since (p\ is a linear functional belonging to Vs 
acting onto an operator \0) belonging to Vq ■ 

In order to see how decoherence works from the new approach, let us consider the simplest case, a quantum system 
whose Hamiltonian has a continuous spectrum: 

H \lü) = lü \lü) lü G [0,oo) (2) 

where lü and oj) are the generalized eigenvalues and eigenvectors of H respectively. The CSCO of this system is just 
{H}. A genèric observable \0) can be expressed in terms of the eigenbasis {\lü)(lü'\} as: 

\0)= í í Ò(LÜ,Lü')\Lü)(Lü'\dLüdLü' = í í Ò{LÜ- : Lü')\Lj;u')dLüdLü' (3) 



where \lü;lü') = \lü)(lü'\ and 0(lü,lü') represents the coordinates of the kerncl \0). The Hamiltonian in the eigenbasis 
{\u;üj')} reads: 

H= lü \lü)(u\ düj = cu S(lü — lü') \lü;lü') dut doj' (4) 



Then, lü5(lü — lü') must be one of the Ò(lü, lü'), since H is one of the observables belonging to Vo ■ Moreover, all the 
observables which commute with H and share the eigenbasis {|cj;w')} must have the following form: 

\0)= J J 0(Lü)\Lü)(Lü\dlü = J J 0(Lü) S(LÜ - Lü') \LÜ-Lü')dLüdlü' (5) 

where 0{lü) supplies the vàlues of the components of \0) in the basis {\lü;lü')}. Therefore, 0(lü)5(lü — lü') must be 
one of the 0(lü,lü'). But, of course, we need also observables which do not commute with H and whose Ò(lü,lü') are 
different than 0(lü) 5(lü — lü'); then, with no loss of physical generality we can postulate as a general case: 

Ò{lü,lü') = 0(lü)6(lü-lü')+0(lü,lü') (6) 

where 0(lü,lü') is a regular function whose precise mathematical properties are listed in g6)i. Therefore, a genèric 
observable \0) reads (see |Ï7|): 

\0)= í 0{lü) \Lü)duj+ í í 0(lü,lü') \ui;ui')du)du)' (7) 



1 Since any singular kernel can be approximated by a regular one, the 0{uj) 8{üj — lü') + 0(lü,u>') are dense en the set of the 
0(üj,u)'). Therefore, we do not loose physical generality, in the sense that the O (ui) 5 (ui — lü') + 0(lü,lü') have all the required 
physical properties up to any order and, then, they are experimentally indistinguishable from the 0(lü,lj'). 
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where \uo) = \u)){w\ and \lo;lo') — are the generalized eigenvectors of the observable \0). We will call the first 

term of the r.h.s of eq.(7) Os (the singular part of the observable \0)), and the second term of the r.h.s of eq.(7) On 
(the regular part of the observable \OY) 

The observables \0) of the form (0) define what we will call "van Hove space", Vq H C V q ; {\uj), \u>;u>')} is the 
basis of Vq H . On the other hand, states are represented by linear functionals belonging to a space Vg , which is 
the dual of Vq H ; therefore, a genèric state (p\ can be expressed: 

(p\ = J p(to) (lo I du + J J p(cj,u') (uo;uo' \duo dio' (8) 

where p(iú,u>') is a regular function, and p(to) and p(co,u)') satisfy the properties p > 0, (p\I) = 1 (where |/) is the 
identity operator) andj-those listed in ]ï^| . {(lo | , (lo; lo' | }, the basis of Vg H , is the cobasis of {\uo), \lo; lo')} defined by 
the following relationstt 

(lo\lo)=S(lo-u>') {lü;lü"\lü';lo"') = 6{lü~lo") S(lj' -lo'") [lo \lo' ; lo" ) = (9) 

Given the expressions @) and for \0) and (p\ respectively, decoherence follows in a straightforward way. Ac- 
cording to the unitary von Neumann equation, the evolution of (p\ is given by: 



(p(t)\ = I p (lo) (lo I dio + I I p(lü,lü') e~ í(lW)í (lo; lo' \ duo duo' (10) 

Therefore, the mean value of the observable \0) in the state (p(t)\ reads: 

(O)p(t) = (p{t)\0) = J p(uo) O(uo) duo + J J p(uo,uo') e-* u - u "> t O(u>,u')du><L·' (11) 

Since p(lo,lo') and 0(lo,lo') are regular functions (see Q for details), if we take the limit for t — > oo, we can apply the 
Riemann-Lebesgue theorem, according to which the second term of the right hand side of the last equation vanishcs. 
Therefore: 



lim (0) p(t) = lim (p(t)\0) = / p(uo) O(uo) duo (12) 

But this integral is equivalent to the mean value of the observable O in a new state (p, | : 

{p*\ = J p(uo) (uo\ duo (13) 

where the off-diagonal terms have vanished. Therefore, we obtain the weak limit: 

lim (0) p(fí = (0) p , (14) 

The next step is to study the formalism undcr the Wigner transformation " symb" . Everything behaves in the usual 
way for the regular parts of \0) and (p\, since these parts satisfy the hypotheses of papers p5[ |. The problem consists 
in defming the Wigner transformation for the singular parts. The singular parts of observables and states read: 



O s = J O(uo) \uo) duo = O(H) p s = I p(uo) (uo \ duo (15) 
Therefore Os is a function of the Hamiltonian: 

H = J uo \lo) duo 

Using the well know properties of the Wigner integral we have that: 



2 These are the generalization of the relations between the basis {\i)} and the cobasis {(j\ } in the discrete case: (j \ i) = Si 
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symbOs = Of(q,p) = 0(H w (q,pj) + 0( — ) (16) 

where S is the characteristic action of the system. Then, in the particular casc whcrc O(uj) = 5(lu — lo') we have from 
cq. ©: 



symb\u'){üj'\ = S(H w (q,p) - u') (17) 
where we have disregarded the 0(^j) as we will always do below. Moreover, for regular functions: 

(p\0) = (symbp\ S ymbO) = J p W \q,p) W \q,p) dq dp (18) 

We will adopt the same equation for the singular parts. This will allow us to define symbps as satisfying: 

(symbp s \symbO s ) = (ps\O s ) (19) 

In doing so we are repeating what we have done at the quantum level when we defined the functional (p|. 
^From eq. ( p7| ) we know that: 

symb\oj') = S(H w (q,p) - uj 1 ) (20) 



It is clear that we cannot normalize this function with the variables and in the domain of integration of (|18j). In 
fact, using the canonical variable —i, the canonical variable conjugated to H, for any function like symb\uj',p') = 
5(H W ($) — ui'), which is a constant for t, the integral will turn out to be infinity. So these functions f(H) are not 
classical densities since they do not belong to L±, if defined integrating over the whole phase space ,and they must be 
normalized in a different way. This fact is not surprising since they are singular functions. But let us observe that, in 
general: 

/•oc />oo 

O s (<j>)=symb 0{lü) \uj) deu = / 0{u) ó(H w (<f>) - u') dui (21) 



which is a function independent of —t, and can be normalized (if necessary) imposing the following conditions: 

i. -We integrate only over the momentum space H (i.e. not over — í), precisely: 

\\O s {<i>)\\= J dH \0(w,p)\S(H -uj')(L· = J \0(H)\dH (22) 

ii. - We chose the regular function O(oj) in the space L\ of the momentum, that is to say: 

\0(u)\dw < oo (23) 



So, we will normalize all the f(H) (if necessary) in this way, and we will perform all the integrations in the Cs space 
by this method. In particular we will use this way of integration when defining functionals of C' s 
Then, in order to satisfy eq.(^), necessarily: 

pZ(q,p) = symb(Lü'\ = 5(H w (q,p) - lo') (24) 

a result already obtained in pi (eqs.(34) and (35)) with different methods. So, PcfuiliP) corresponds to a density 
function extremely peaked over the classical trajectory defined by the conservation law H w (q,p) = lo'. 
^,From eq. (|Ï5|) : 

/•oo />oo 

pj(q,p)= P (u)pZ(q,p)cL·= p(u>) S(H W (q,p) - u>') du = p(H w (q,p)) (25) 
Jo Jo 

and it is a constant of the motíon, as it was expected. Moreover, p ] g(q,p) > since ps(^) > 0. This means that 
pg (q,p) is the statistical ensemble of the density functions p^il'P) extremely peaked over the classical trajectòries 
defined by the conservation law H w (q,p) — cu' and weighted by the probabilities p(lo). From eq. (|14|) we know that 
only the singular part must be considered after decoherence. Thus, for t — > oo the classical density pg(q,p) is resolved 
as a set of classical trajectòries. 
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At this point, we have defined the Wigner transformation both for the regular parts and for the singular parts of 
observables and states. As the result, the singular parts share the same usual properties with the regular parts, since 
we have postulated such properties to define the Wigner transformation of the singular parts. In fact, these usual 
properties follow from eq.(16) (see [Ï8| for details). E.g., eq.(12) is also vàlid when the Wigner transformations are 
involved: 

Hm (p w \O w ) = [ pf{H) Oj(H) dH (26) 

í^oo J 

This means that any observable \0) becomes OÍj 7 (q,p) in the limit t — ► oo, and behaves in a classical way. 

After this presentation of the formalism in the simplest case, let us study the general case. In general, we must 
consider a CSCO, {H, Oi, O n }, whose eigenvectors are \oj, 01, o n ). In this case, will be diagonal in the 
variables us, ui' but not in general in the remaining variables. Therefore, a further diagonalization of is necessary: 
as the result, a new set of eigenvectors {\ui,pi, ...,p n }}, corresponding to a new CSCO {H, Pi, P n } emerges. This 
set defines the eigenbasis ...,p n ), \w,pi, ...,p n ;ui' ,p[, ...,p' n )} of the van Hove space of observables Vq H , where: 

\u,Pi, ...,p n ) = \w,Pi, -,Pn) (u,Pi, -,Pn\ (27) 

\UJ,P1, ...,p n ;üj',p[, ...,p' n ) = \v,Pl, ...,Pn) {^',p'l, -,Pn\ 

(p*| will be completely diagonal in the cobasis of states {(cc>,pi, ...,p n \, (u>,Pi, ■■-,p n \ lo' ,p' lli: .,p' n \} corresponding to 
the new eigenbasis of Vq H (see |ï^] for details). And, most important, for this system eq.(|25|) reads: 

/■OO />00 

Ps(l,P)=J2 P(")pLp 1 ,-, P Jl>P) d " = Y, Pn6(H w (q,p)-u J i)5(Pr(q,p)-p>)du, (28) 
p Jo p Jo 

where the density p^ pi Pn corresponds to a density function extremely peaked over the classical trajectory defined 
by the conservation laws: 

H w (q,p)=cü', pW(q,p)=pi, ... P™(q,p)=p n (29) 

Therefore, for t — > oo the classical density p^f (q,p) is resolved again as a set of classical trajectòries. 

As this presentation shows, decoherence does not require the interaction of the system of interest with the environ- 
ment: a single closed quantum system can decohere. The diagonalization of the density operator does not depend on 
the openness of the system but on the continuous spectrum of the system's Hamiltonian. This means that the prob- 
lem of providing a general criterion for discriminating between system and environment vanishes in the self-induced 
approach. This fact leads to an additional advantage of the new way of conceiving decoherence. As we have seen, 
in many cases the einselection approach requires to introduce assumptions about the observables which will behave 
classically in order to decide where to place the boundary between system and environment. The new approach, 
on the contrary, provides a mathematically precise definition of the observables regarding to which the system will 
decohere. In fact, there are two kinds of such observables: 

a) Observables that commute with the Hamiltonian, which are represented by the singular kernels 0(u>) S(lu — lo'). 

b) Observables that do not commute with the Hamiltonian, which are represented by the kernels 0(u>) S(lü — lo') + 
0(w,u>'), where 0(u>,u>') is a regular function. In other words, these observables have a regular part 0(u>,u>') and a 
singular part 0(u>) S(lu — uj') in the eigenbasis defined by the system's Hamiltonian. 

This definition is completely general and does not require to introduce any prior assumption about the classical 
behavior of certain observables. 

When the phenomenon of decoherence is viewed from this new perspective, it does not need to be conceived as "a 
justification for the persistent impression of reality" |ÏJ . Classicality is not a perceptual result of the correlations 
between the observed system and the observer's brain though the environment: the emergence of classicality is a 
consequence of the own dynamics of a closed quantum system. In other words, from the self-induced approach, 
decoherence is a relevant element for explaining the emergence of classicality, not our perception of classicality. 

IV. DECOHERENCE IN A CLOSED UNIVERSE 

If the transition from quantum to classical does not require the split of the universe into subsystems as a necessary 
condition, then decoherence can take part in the account of how the universe as a whole behaves classically. In 
this section we will apply the self-induced approach to a simple quantum-cosmological model in order to show how 
classicality arises in this case. 
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A. The model 



Let us consider the flat Roberson- Walker universe ( |ïï| , [^0| ) with a mètric: 

ds 2 = a 2 (r)) (di] 2 - dx 2 - dy 2 - dz 2 ) (30) 

where 77 is the conformal time and a the scale of the universe. Let us consider a free neutral scalar field $ and let us 
couple this field with the mètric, with a conformal coupling (£ = ~). The total action reads S — S g + S / +Si, and 
the gravitational action is: 

S g = M 2 J dn[-\a -V(a)] (31) 

where M is the Planck mass, a= da/di], and the potential V contains the cosmological constant term and, eventually, 
the contribution of some form of classical matter. We suppose that V has a bounded support < a < a±. We expand 
the field <£> as: 

/+00 
/ k (r?)e- íkx dk (32) 
-00 

where the components of k GR 3 are three continuous variables. 
The Wheeler-De Witt equation for this model reads: 

i/*(a,$) = (h g + h f + h t )^{a,<S>) = (33) 

where: 

h 9 = ^d 2 a + M 2 V(a) 



\J{dl~k 2 f 2 ) dk 



h i = ^m 2 a 2 í fidk (34) 

with m the mass of the scalar field, k/a the linear momentum of the field, and —d/df^. 
We can now go to the semiclassical regime using the WKB method ( |2Ï|), writing ^(a, $) as: 

*(o, $) = exp[iM 2 S*(a)] X (a, $) (35) 

and expanding S" and x as: 

5' = S' + M- 1 5 1 + ..., x = Xo + A/- 1 xi + ... (36) 



To satisfy eq. (33) at the order M , the principal Jacobi function S'(a) must satisfy the Hamilton- Jacobi equation: 

2 



dS_ 

da 



2V(a) (37) 



We can now define the (semi) classical time as a parameter 77 = T)(a) such that: 

d dS d 1——-. . d ,„„ N 

j- = "T J" = ±V2^(a) — (38) 
dr\ da da da 

The solution of this equation is a = ±F(rj,C), where C is an arbitrary integration constant. Different vàlues of this 
constant and of the ± sign give different classical solutions for the geometry. 

Then, in the next order of the WKB expansion, \ satisfy a Schròdinger equation that reads: 
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where: 



precisely: 



where: 



1 % = h(jí)X 



h{rj) = h f + hi(a) 



m = - 2 



d 2 

-qJ2 +^k(a)/k 



r/k 



0|(a) = íl^(a) = m 2 a 2 + A: 2 = m 2 a 2 + w 



(39) 



(40) 



(41) 



(42) 



where tu — k 2 and k = |k|. So the time dependence of the Hamiltonian comes from the function a = a{rj). 

Let us now consider a scale of the universe such that a ou t 3> a±. In this region the geometry is almost constant. 
Therefore, we have an adiabatic final vacuum |0) and adiabatic creation and annihilation operators and a^. Then 
h = h(a ou t) reads: 



h = J il^a^a^dk 

We can now consider the Fock space and a basis of vectors: 

|ki,k 2) ..,k n )Sí|{A}) =aí 1 4 a ...4 n ...|0> 



(43) 



(44) 



where we have called {k} the set ki, k2, k„, where eventually n goes to infinity. The vectors of this basis are 
eigenvectors of h: 



h\{k}) =Lü\{k}) 



where: 



w= (m 2 a 2 out +zu)^ 
ke{k} ke{k} 



We can now use this energy to label the eigenvector as: 

IW)=K[k]) 



(45) 
(46) 

(47) 



where [k] is the remaining set of labels necessary to define the vector unambiguously. {\u>, [k])} is obviously an 
orthonormal basis, so eq. (|43]) reads: 



u\u,\k])(u, \k]\dud\k] 



(48) 



B. Decoherence in energy 

In this case, a genèric observable \0) £ Vq H reads (see eq.(^)): 

\0) = J O(o;,[k],[k']) |w,[k];[k'])dwd[k]d[k / ] + J J 0(ui, [k],u/, [k']) \u, [k];u\ [k'])duid[k] dui' d[k'] 

where: 

k [k];[k']) = k [k])k, [k']|) | W , [k]; J, [k']) = | Wj [k])<ü/, [k']| 
and a genèric state (p\ G can be expressed as (see eq.(||)): 



(49) 
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(p\ = J p(w,[k] ( M)( W) [k];M| <L·d[k]d[k'] + I J p(u,[kW,[k']) (cj í [k] ] oj',[k']\du,d[k]cL·'d[k'] (50) 

where {(w, [k];[k'] | , [w, [k]; w' , [k']|}, the basis of ^ s yíí , is the cobasis of {\lü, [k];[k']), \u, [k];o/, [k'])}. Thcn, the mean 
value of the observable |0) in the state (p(í)| reads: 



(0) p(t) = (p(t)\0) = J p( W ,[k];[k'])0(^[k];[k'])^d[k]d[k']+ (51) 

+ J J p(uj, [k],uj', [k']) e -*(<"-«0* 0(w, [k],c/, [k']) du d[k] dw' d[k'] (52) 
Taking the limit for í — > oo and applying the Riemann-Lebesgue theorem, we obtain: 

lim (0) m = hm (p(t)|0) = / p( W ,[k];[k'])0( W ,[k];[k'])^d[k]d[k'] (53) 

i— ► oo v ' t—*-oo J 

And this integral is equivalent to the mean value of the observable \0) in a new state (p*|: 

(p*| = / P(w, [k];[k']) (w, [k];[k'] | dwd[k] d[k'] (54) 

This new state (p*| is the equilibrium time-asymptotic state, which is diagonal the variables u>,lj' as decoherence in 
energy requires. 

C. Decoherence in the remaining dynamical variables 

In this case, (p*| is diagonal in the variables uj,uj' but not in the remaining variables. This means that a further 
diagonalization is necessary. 

Let us observe that, if we use polar coordinates for k, eq.(|32"|) reads: 

$(x, n) = J falm dk (55) 

Irn 

where: 

<pkim = fk,i(r),r)Yl(6,<p) (56) 

where fe is a continuous variable, Z = 0,1, ; m = —l, 1; and are spherical harmònic hmctions. So the indices 
fe, l, m contained in the symbol k are partially discrete and partially continuous. 

As (pt\ = (p*|, then p*(u), [k],[k']) = p(ui, [k],[k']) and, therefore, there exists a set of vectors {\oj, [1])} such that: 

J p{u, [k],[k>]) | W) d[k>] = p(iu, [1]) \u>, [l]) [k] (57) 

namely, {\ui, [1])} is the eigenbasis of the operator p(ui, [k],[k']). Then p(ui, [1]) can be considered as an ,e>rdinary 
diagonal matrix in the discrete indices / and m, and a generalized diagonal matrix in the continuous index fets. Under 
the diagonalization process, eq.(|54|) is written as: 

(P*\ = I U$ [!] = [!']) <í U(ff (lü, [l"}, [1"'] | U^' ] cL·dlk] d[k>] d[l] d[l>] d[l"] d[l'"\ (58) 

where uffl is the unitary matrix used to perform the diagonalization andEl: 



s E.g.: We can deal with this generalized matrix by rigging the space 



this procedure allows us to define a generalized state eigenbasis. 
diagonalize Hamiltonians with continuous spectra in [^3), |MJ . 



V/ H and using the Gel'fand-Maurin theorem 1 

It has been usec 



br systems with continuous spectrum. 
5|, etc. 



<5([1]— [!']) is a Dirac delta for the continuous indices and a Kronecker delta for the discrete ones. 
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p(u,[lUl'})=p(u>,[l})8([ï\-[l'}) (59) 

Sincc: 

p(u, [1],[1]) = p(u, [1]) = I U$ p(u, [k],[k']) t/Jj"' 1 d[k] d[k'} (60) 

we can define: 

(u, [1]| = (c, [1], [1]| = | C/M ( W , [k], [k']| üjpt dpc] d[k'] (61) 

We can repeat the procedure with vectors (u>,lo', [k], [k']| and obtain vectors (u>,lo', In this way we obtain a 
diagonalized cobasis {(ui, [1]|, (u>,lu', [1]|}. So we can now write the equilibrium state as: 

(p,| = yp(a;,[l])(w,[l]|^d[l] (62) 

Since vectors (w, [1]| can be considered as diagonal in all the variables, we have obtained decoherence in all the 
dynamical variables. 



D. Emergence of trajectòries 

Let us restore the notation {1} = (ui, [1]), {k} — (ui, [k]) as in eq.([44]), and let us considcr the coníiguration kets 
\{x}) — [x]}. Since we are considering the period when a ~ a ou t, the system with Hamiltonian ( [43] ) is just a set of 
infinite oscillators with constants Qk(a ou t) that represents a scalar field with mass ma ou t. Then, we are just dealing 
with a classical set of N partides, with coordinates [x] and momenta [k]. Therefore, we can introduce the Wigner 
function corresponding to the generalized state \{l})'- 

^([x], [k]) = n- iN I ({í}|x + A)(x - A|) e *WM é n \ (63) 

Using the same reasoning that we have used to obtain eq. (p8|) (see also paper p8| ) p|^([x], [k]) reads: 

p™([x},[k]) = l[5(L™([x},[k])-l i ) (64) 

i 

where ([x], [k]) is the classical observable obtained from Li (that corresponds to indices 1) via the Wigner integral 
(considering h = Lq and including among the indices i). Now, with the new notation eq.(|62j) reads: 

ÍP*\ = I P{i}({l}\d{l} (65) 

Then, if we call: 

pf ([x], [k]) = ir- iN J (p,|x + A)(x - A|) e^W'M d in \ (66) 

we obtain: 

pf([x], [k]) = pf(C([x], [k]), Lf([x], [k]), ...) (67) 

So finally: 

pf([x],[k])~ | m pr([x],[k])<5({L^}-{0)d{0 = / p m \X[6(LY -k)d{l] (68) 

The last equation can be interpreted as follows: 

i.- 6({L W } — {l}) is a classical density function, strongly peaked at certain vàlues of the constants of motion {/}, 
corresponding to a set of trajectòries, where the momenta are equal to the eigenvalues k (i = 0, 1, 2, ...). 
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ii. - p{;} is the probability to be in one of these sets of trajectòries labelled by {l}. Precisely: if some initial density 
matrix is given, from eq.(p5|) it is evident that its diagonal terms p^ are the probabilities to find the density function 
6({L W } — {l}) in the corresponding classical equilibrium density function pj/ ([x], [k]), namely, the probabilities of 
the trajectòries labelled by {1} = (w, [1]). n 

iii. - Let a be the coordinate classically conjugated to 1 and let be ao the coordinate a at time rj = 00, then we obtain 
the classical trajectòries: 

a i= U^+a oi (69) 

iv. - Let us now call p^j — piau^]- Actually, P{/}[ ao ] is not a function of ao; it is simply a constant in ao, since ao 
is only an arbitrary point and our model is spatially homogenous. Then we can write: 



P{1}[, 



/7l 
P{i}l*o] n s i a i - a 0l )d[a ] (70) 



In this way we have changed the rolc of ao: it was a fixed (but arbitrary) point, and now it is a variable that moves 
all over the space. Then eq.(^) reads: 

r, n n 

p^Qx], [k]) ~ j p W[ao] l[6(Lf - k) J] S( aj ao,) d[ao] d{l} (71) 



So, if we call 



^ }[ao] ax], pe]) = n<w - u) n ^ - a ^-) ( ?2 ) 

í=0 3=1 



we have: 

«W/m run / m 



^(MJkj) ~ J p {I}[aa] ^ }[ao] ([x],[k])d[ao]d{l} (73) 

^From eq.([7^) we see that p|^r a ,([x], M) ^ om Y m a narrow strip around the classical trajectory ( p9| ) defined by 
the momenta {1} and passing through the point [ao] (actually the density function is as peaked as it is allowed by the 
uncertainty principle; its width is essentially a 0(4), since the S— functions of all the equations.. Decome Dirac's deltas 
only when % — ► 0). Therefore, eq.(|73]) describes the classical behavior of our model of universeE. 

Let us sum up the main steps of our argument. When r\ — > oo, the quantum density p becomes a diagonal density 
matrix p*. The corresponding classical distribution p^ ([x], [k]) can be expanded as a sum of classical trajectòries 
density functions P|j}[ ao ]([x], [k]), each one weighted by its corresponding probability P{i}[ ao ]. Going back to eq. ([73|) , 
prnr ao j([x], [k]) = la i n }[ ao ]([ x ], H) is the density corresponding to the set of n points (let us say, galàxies), each 



one of them moving over a trajectory defined by eq.(69), where eventually n goes to infinity. So, as the limit of our 
quantum model we have obtained a classical statistical mechanical model, and the classical realm appears. 

In summary, we have proved that the density operator is translated into a classical density, via a Wigner function, 
and it is decomposed as a sum of densities peaked around all possible classical trajectòries, each one of these densities 
weighted by their own probability. Therefore, our quantum density operator behaves in its classical limit as a statistical 
distribution among a set of classical trajectòries. Similar results are obtained in papers pq] and ETfl. 



5 We could also add a non-relevant equation, something like t = íüt] + to, which would define a choice of our clock's time. 

6 In this section, as in section IIB, we have faced the following problem: pj^([x], [k]) is a a constant that we want to decompose 
in functions P{n[a ] (M' M) w hich are difFerent from zero only around the trajectory (^) and therefore are variables in a. Then, 
essentially we use the fact that if f{x,y) — g(y) is a constant function in x, we can decompose it as: 

g(y) = J g(y)S{x - x )dx 

namely, the densities S(x — xo) are peaked in the trajectòries x = xo — const., y = var. and, therefore, are functions of x. This 
trajectòries play the role of those of eq.(|70[). 

As all the physics, including thecorrelations, is already contained in eq. (p^) , the reader may just consider the final part of 
this section, from eq.djd) to eq.(|73[) a didactical presentation. 
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V. OPERATORS AND FIELDS 



Up to this point, all our reasoning was made in the context of the Schròdinger picture, as it is usual in quantum 
mechanics. But in quantum field theory, on the contrary, the Heisenberg picture is the usual scenario. However, it is 
quite easy to reformulate the argument in the new picture, since our main equations are (14) and (26) for the general 
case, and (53) and its corresponding classical version for the cosmological case. In both cases only the singular parts 
Os and ps = p* appear, and both functions are time-independent. Then, whereas in the Schròdinger picture we 
have that O = O(ío) is time-constant and p(t) — > = ps is time-constant in the limit, in the Heisenberg picture 
p = p{to) is time-constant and 0(t) — > O* — Os is time-constant in the limit; nevertheless, both results coincide 
since, when t — » oo, = ps and = Os- Then, if wc consider the Heisenberg picture, when t — > oo the operators 
/k(?y) of eq.(32) become singular and also equal to the constants fu; therefore, the corresponding are constant 
and classical. The same can be said for the field $(77, x). of eq.(32), that becomes a time constant $(x). 

With this procedure, in certain sense we have lost the dynamics of the field. However, it can be recovered if we 
compute a "master equation" from any asymptotic expansion of the Riemann-Lebesgue theorem or if we make an 
analytic continuation in the Liouville complex energy plane as in paper [p8| . Such a master equation coincide with 
the usual one and also with the Lindblad approach. 



VI. CONCLUSION 



In this paper our aim was to rigorously present the self-induccd approach to decohcrence, according to which the 
phenomenon of decoherence is the result of the own dynamics of a closed system governed by a Hamiltonian with 
continuous spectrum. From this approach, the interaction between the system and the environment is not needed and, 
therefore, a single closed quantum system can decohere. This feature makes this approach particularly appropriate 
for addressing problems in quantum cosmology, since the universe is, by definition, a quantum closed system with no 
environment to interact with. In this context, we have applicd the self-induced approach to a quantum cosmological 
model, showing that decoherence in energy and in the remaining dynamical variables obtains with no reference to an 
environment and without assuming in advance which observables will behave classically. 
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